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ABSTRACT 

We analyze the scaling of the X-ray power density spectra with the mass of the 
black hole on the example of Cyg X-1 and Seyfert 1 galaxy NGC 5548. We show that 
the high frequency tail of the power density spectrum can be successfully used for 
determination of the black hole mass. We determine the masses of the black holes in 6 
Broad Line Seyfert 1 galaxies, 5 Narrow Line Seyfert 1 galaxies and two QSOs using 
available power density spectra. The proposed scaling is clearly appropriate for other 
Seyfert galaxies and QSOs. In all but 1 normal Seyferts the resulting luminosity to the 
Eddington luminosity ratio is smaller than 0.15, with a source MCG -6-15-30 being 
an exception. The applicability of the same scaling to Narrow Line Seyfert 1 is less 
clear and there may be a systematic shift between the power spectra of NLSl and SI 
galaxies of the same mass, leading to underestimation of the black hole mass. However, 
both the method based on variability and the method based on spectral fitting show 
that those galaxies have relatively low masses and high luminosity to the Eddington 
luminosity ratio, supporting the view of those objects as analogs of galactic sources in 
their high/soft or very high state based on the overall spectral shape. Bulge masses 
of their host galaxies are similar to normal Seyfert galaxies so they do not follow 
the black hole mass - bulge mass relation for Seyfert galaxies, being evolutionary less 
advanced, as suggested by Mathur (2000). The bulge mass-black hole mass relation 
in our sample is consistent with being linear, with black hole to bulge ratio ~ 0.03 % , 
similar to Wandel (1999) and Laor (1998, 2001) for low mass objects but significantly 
shifted from the relation of Magorrian et al. (1998) and McLure & Dunlop (2000). 

Key vi^ords: galaxies: active - accretion, accretion discs - black hole physics - binaries 
- X-rays: stars - galaxies: Seyfert, quasars - X-rays. 



1 INTRODUCTION 

Determination of the masses of black holes in active galactic 
nuclei (AGN) is a key element in studies of the nature of 
accretion process and evolution of an active nucleus. 

Several methods were used so far for this purpose: esti- 
mation of the virial mass from the emission line data through 
reverberation techniques (see e.g Wandel 1999), estimation 
of the virial mass on the basis of the source luminosity and 
estimates of the ionization parameter in Broad Line Region 
(e.g. Wandel, Peterson & Malkan 1999), and model fits (usu- 
ally of accretion discs) to the broad band spectra (e.g. Edel- 
son & Malkan 1986, Sun & Malkan 1989, Merloni, Fabian & 
Ross 2000). 

Suggestion that variability properties of the sources are 
affected by their mass and/or luminosity was made by Barr 
& Mushotzky (1986) and Wandel & Mushotzky (1986), and 



subsequently explored by many authors (e.g. McHardy 1989, 
Green, McHardy & Lehto 1993, Turner et al. 1999, Leighly 
1999a). 

A convenient and powerful method of determination 
of black hole mass from variability properties was recently 
used by Hayashida et al. (1998) and Hayashida (2000). The 
method uses the normalization of the X-ray power density 
spectra (hereafter PDS) in comparison with Cyg X-1 in the 
hard state. The method, however, lead to surprisingly small 
values of the black hole masses and, consequently, to lumi- 
nosity to the Eddington luminosity ratios frequently larger 
than one. 

In the present work we reconsider the method of 
Hayashida et al. (1998). We determine better scaling co- 
efficient on the basis of the new available power spectra of 
Cyg X-1, the new power spectrum of NGC 5548, and the 
independent black hole mass measurement for this AGN. 
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We determined the black hole masses for a number 
of AGN with the PDS available in the literature. We in- 
cluded both the classical Seyfert 1 galaxies (i.e. Broad Line 
Seyfert 1 galaxies; hereafter BLSl) and Narrow Line Seyfert 
1 galaxies (NLSl), characterized by much narrower broad 
lines {Vfwhm < 2000 km s"\ Osterbrock & Pogge 1985). 

We discuss the consequences of those black hole mass 
determination for the accretion rate in BLSl and NLSl 
galaxies and for their evolutionary status. 



2 METHOD 

2.1 Determination of black hole mass in NGC 
5548 

Seyfert 1 galaxy NGC 5548 (z = 0.0174) is one of the most 
frequently studied AGN. The mass of the central black hole 
in this object is therefore determined relatively accurately. 
Peterson & Wandel (1999) give the value G.9,t\l x W Mq 
from reverberation studies of the various lines emitted by 
the Broad Line Region. This technique may potentially suf- 
fer from systematic errors up to a factor of 3 (Krolik 2000) 
but recent comparison of statistical properties of reverber- 
ation and dynamical mass measurement indicated smaller 
errors (Gebhardt et al. 2000). Virial mass of the black hole 
in NGC 5548, determined on the basis of ionization method 
for a number of emission lines is equal to 5.9 ± 2.5 x IQ^ Mq 
(Peterson & Wandel 2000). 

Those measurements confirmed earlier results based 
on the study of the BLR which indicated the mass about 
lO^M© (Wanders et al. 1995), or somewhat smaller (Done 
& Krolik 1996). Rokaki et al. (1993) argued for the value 
between 5 x lO^Af© and 6 x lO^'M© and Kaspi et al. (2000) 
obtained values 9.4 - 12.3 x lO^'M©. 

Fits of a simple disc/corona model to the optical/UV 
continuum data from the AGN Watch (Clavel et al. 1991, 
Peterson et al. 1992) favored the value 10* M© although the 
value 6x lO^M© was also acceptable (Loska & Czerny 1997). 
Fits of the more advanced disc/corona model of Witt, Cz- 
erny& Zycki (1997) which was based also on the X-ray data 
indicated rather larger value 1.4 x 10* M© (Kuraszkiewicz, 
Loska & Czerny 1997). 

The value of the black hole mass 6.8^10 x lO^M© seems 
to be determined most accurately, since it resulted from the 
analysis of several emission lines coming from the Broad Line 
Region. It is also consistent with other determinations, tak- 
ing into account quite large and unspecified errors in mea- 
surements based on other methods. 



2.2 Power density spectrum of NGC 5548 and the 
comparison Avith galactic sources 

Time variability of the X-ray and optical emission from the 
nucleus of NGC 5548 was studied by a number of authors 
(e.g. Papadakis & Lawrence 1993, Clavel et al. 1992, Czerny, 
Schwarzenberg-Czerny & Loska 1999, Chiang et al. 2000 and 
references therein). The power density spectrum is basically 
featureless, as in other AGN, and the variability is caused 
by some stochastic process (Lawrence et al. 1987, McHardy 
& Czerny 1987, Czerny & Lehto 1997). Since Comptoniza- 
tion is directly responsible for the formation of the hard 
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Figure 1. Power density spectrum times frequency for NGC 5548 
from Chiang et al. (2000) (stars connected with continuous line) 
and for galactic black holes: Cyg X-1 in transition state from Cui 
et al. (1997) (open triangles connected with continuous line), Cyg 
X-1 in hard state from Lin et al. (2000) (solid triangles connected 
with long dash line), GRS 1758-258 from Lin et al. (open octagons 
connected with short dash-long-dash line), IE 1740.7-2942 from 
Lin et al. (stars connected with dot-short-dash line), GX 339- 
4 from Lin et al. (solid squares connected with dot line). Power 
density spectrum of NGC 5548 was shifted by 6.83 in log{y) which 
corresponds to rescaling this source to a black hole mass of IOM0 . 



X-ray spectra this variability must be related to variable 
seed photon fiux and/or hot plasma properties but the exact 
physical mechanism is still unknown (for a review of X-ray 
spectra and variability of AGN, see e.g. Mushotzky, Done & 
Pounds 1993, Nandra et al. 1997, Leighly 1999a,b; but see 
also Abrassart & Czerny 2000 and the references therein for 
an alternative view). The same mechanism is expected to 
operate in galactic black holes (GBH; for a review, see e.g. 
van der Klis 1995, Cui 1999, Poutanen 2000). 

The basic similarity of AGN and GBH with respect to 
the spectra and variability was discussed by a number of 
authors (e.g. Zdziarski 1999). There are, however, also some 
systematic differences refiecting the direct dependence on 
the mass of the central body: (i) the temperature of the 
disc component is about two orders of magnitude higher in 
galactic sources than in AGN (ii) GBH in their soft/high 
state are generally less variable than AGN dominated by 
the disc component. 

Recent study of NGC 5548 with RXTE satellite allowed 
to determine an exceptionally accurate PDS for this source 
in the X-ray band (Chiang et al. 2000). 

We compare PDS of NGC 5548 with PDS of a number 
of galactic sources in various luminosity states (see Fig. m. 
We plot these spectra in the form of z^ x power for more con- 
venient comparison. We rescale the PDS of NGC 5548 to the 
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Figure 2. Power density spectrum times frequency for NGC 5548 
from Chiang et al. (2000) (histogram: continuous line - XTE, 
dotted line - EUVE; stars connected with continuous line - fit) and 
for galactic black holes: Cyg X-1 in hard state from Revnivtsev 
et al. (2000) (open squares connected with dot line), Cyg X-1 in 
soft state from Gilfanov et al. (solid triangles connected with long 
dash line). Power density spectrum of NGC 5548 was shifted by 
6.83 in logiy) which corresponds to rescaling it to a black hole 
mass of IOMq. 



mass of the black hole IOMq, appropriate for Cyg X-1 (see 
Nowak et al. 1999 and the references therein) through simple 
horizontal shift by a factor of 6.8 x 10^. We show separately 
the comparison of NGC 5548 with the high quality data of 
Cyg X-1 in its high/soft (Gilfanov et al. 2000) and low/hard 
(Revnivtsev, Gilfanov & Churazov 2000) states (see Fig. m). 

The exact comparison of the two distributions is how- 
ever difficult since the quality of the PDS for NGC 5548 is 
much lower than that of Cyg X-1. The fit to the overall shape 
of the PDS found by Chiang et al. (2000) suggests that there 
is a systematic shift between the two distributions (see con- 
tinuous line in Fig. Q) by about 0.4 in the logarithmic scale. 
The high frequency turn off is best reproduced by Cyg X- 
1 in soft or in the transition state, but the normalization 
is considerably higher, confirming much stronger variability 
than in the soft state of Cyg X-1. On the other hand the 
fit is strongly influenced by the results from EUVE which 
represents different energy band. XTE results are roughly 
consistent with the distribution of Cyg X-1 in the high fre- 
quency part. 

The exact position of the high frequency part of the 
spectrum is the basis of the mass measurement proposed by 
Hayashida et al. (1998). The dependence on the luminosity 
is constrained to low frequencies (Belloni & Hasinger 1990) 
so the high frequency part (above 10 Hz for galactic sources 
and correspondingly lower for AGN) is promising for mass 
determination. 

According to Hayashida et al., the frequency i^o.ooi 



where the power spectrum in i/ x Power representation is 
equal 10""^ scales with the mass, log(Ms_ff/MQ) = C — 
logfo.ooi- Hayashida et al. (1998) adopted the value C equal 
2.66 at the basis of their Cyg X-1 spectrum. The new hard 
state data indicates the value of i^o.ooi by more than a factor 
of 2 higher, (C = 3.1). 

The power spectrum of NGC 5548, if well represented 
by the fit, is slightly shifted towards higher frequencies and if 
the mass of Peterson & Wandel (1999) is adopted, C = 3.5. 

The difference between the value of C from Cyg X-1 
and NGC 5548 is larger than the formal error given on the 
mass measurement of NGC 5548 (0.1; Peterson & Wandel 
1999). 

This difference may be partially caused by the poor 
quality of NGC 5548 data. Simulations indicate that the de- 
termination of the high frequency slope of the PDS may be 
biased by the power leak to high frequencies connected with 
the window function if the PDS slope is close to -2 (Green, 
McHardy & Lehto 1993; see also Czerny et al. 1999). 

The shift may also reflect the uncertainty in our knowl- 
edge of the masses: neither the mass of Cyg X-1 nor of NGC 
5548 is known very accurately, if systematic errors connected 
with reverberation method are taken into account (Krolik 
2000). From this point of view, it would be more convenient 
to use the galactic source GRO J1655-40 because the mass 
determination for this source is by far the most accurate 
(A/ = 7.2 ± 0.22Mq, Orosz & Bailyn 1997; see the discus- 
sion of Ziolkowski 2001). Unfortunately, the currently avail- 
able power spectrum for this source (Sunyaev & Revnivtsev 
2000) is not as accurate as for Cyg X-1 so at the basis of 
this source the value of C should be roughly between 2.9 
and 3.3, supporting the result from Cyg X-1. 

The systematic shift between NGC 5548 and Cyg X-1 
comes predominantly from the EUVE results biasing the fit 
proposed by Chiang et al.(2000). The XTE results them- 
selves (a histogram in Fig. 0) are roughly consistent with 
Cyg X-1. 

Therefore, assuming the value C = 3.1 for both galactic 
objects and AGN seems reasonable. 



3 DETERMINATION OF THE CENTRAL 
BLACK HOLE MASSES IN AGN 

We use the available PDS for a number of Seyfert galaxies in 
order to determine their masses; PKS 2155-304, NGC 4051, 
NGC 4151 and NGC 5506 is taken from Hayashida et al. 
1998, MCG -6-30-15 is from Nowak & Chiang (2000), and 
NGC 3516 from Edelson & Nandra (1999). PDS for NGC 
7469 is from Nandra & Papadakis (2001), but with the nor- 
malization different by a factor of 2 from the original nor- 
malization which assures that the integral of the PDS gives 
the variance. All distributions, together with NGC 5548, are 
plotted in Fig. ||. 

The PDS of various AGN are essentially similar to each 
other and to the PDS of NGC 5548. Horizontal shifts be- 
tween them can be explained as the differences in their black 
hole masses. However, there are also some differences in to- 
tal normalization. Some of the effect may be connected with 
the finite length of the data, not extending beyond the max- 
imum of the PDS at the longest timescales. However, the 
effect may be also connected with their luminosity. Similar 
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4 BULGE MASSES OF HOST GALAXIES 
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Figure 3. Power density spectrum times frequency for NGC 
5548 (stars connected with continuous line), PKS 2155-304 (open 
hexagons connected with continuous line), NGC 3516 (solid 
squares connected with dot line), NGC 4051 (open squares con- 
nected with short dash line), NGC 4151 (solid triangles connected 
with dot-short-dash line), NGC 5506 (open triangles connected 
with long dash line), NGC 7469 (solid hexagons connected with 
short dash line), MCG -6-30-15 (open octagons connected with 
short-dash-long-dash line). PDS of these objects are used from 
papers cite in text (see Section 3). 



dispersion is also seen among the galactic sources. It clearly 
limits the accuracy of the method. 

We determine the black hole masses from their PDS 
using the method of scaling adopted by Hayashida et al. 
(1998). We adopt the PDS value of 10-^(rms)2 Hz-^Hz 
as the characteristic value and we compare the frequencies 
J'o.ooi where these value is reached in various sourc es (t hick 
horizontal line in Fig. w. Our formula is (see Sect. 2.2) 



log(Mflif/MQ) = 3.1 - log!/o.ooi- 



(1) 



For those sources analyzed by Hayashida et al. (1998) 
and Hayashida (2000) for which no better PDS are cur- 
rently available we simply adopt their results after scaling 
the masses by a factor 2.8 (shift by 0.44 in logarithm), which 
results from systematic difference betw een new and previ- 
ous results for Cyg X-1 (see Sect. 2.2). Those sources are: 



PHL 1092, PG 1244, IRAS 13224, IH 0707 from Hayashida 
(2000), and 3C 273, ESO 103-G35, PKS 2155-304, NGC 
4051, NGC 4151 and NGC 5506 from Hayashida et al. 
(1998). 

All masses of central black holes, are given in Table hi. 



We determine the bulge masses of the considered AGN us- 
ing the data available from the literature and following the 
approach of Laor (1998) and Wandel (1999). 

No data are available on the host galaxy of PHL 1092 
and PKS 2155-304, so we were unable to determine bulge 
masses of host galaxies in the case of these two objects. 

For NGC 5548, NGC 4051, NGC 3516, NGC 4151, 
NGC 5506 we take the blue bulge magnitude A'Ib from 
Whittle (1992), we calculate the visual magnitude as 
Mv = Mb - 0.95 (Worthey 1994), bulge luminos- 
ity as [og{Li,uige/Le) = OA{-Mv + 4.83) and finally 
log{Miuige/MQ) = 1.18log{Lb^ige/LQ) - 1.11, after Magor- 
rian et al. (1998). Our values of masses for 5548, NGC 4051, 
and NGC 4151 differ by a factor of ~ 2 from those of Wan- 
del (1999) due to the A'lv = Mb — 0.95 correction applied 
by us. 

For 3C 272 we take the visual magnitude of the host 
galaxy M^?""' = -22.58 from Laor (1998), and the type of 
the galaxy E4 after Bahcall et al. (1997), we then calcu- 
late the bulge luminosity from the formula M^°*"' - My = 
0.324r - 0.054r^ -I- 0.00047t^, where r = T -(- 5 and T is 
the Hubble type of the galaxy and we further proceed as 
described above. 

In the case of remaining AGN, apart from IRAS 13224- 
3809, we took their apparent magnitudes and redshifts from 
NED, we corrected for the contamination of the total flux by 
the nuclear emission using the measured equivalent width of 
the Hg line, EWHp, as M^b"'"' = M^'^''""'''^'' - 2.51og(l - 
1.10 * {1 + z)EWHi3/100), after Whittle (1992), assuming 
the spectral slope equal -1 in all objects. EWHp values 
were taken from Morris & Ward (1988) for MCG -6-30-15, 
and ESO 103-G35, and from Miller et al. (1992) for PG 
1244-f 226, and from Giannuzzo & Stirpe (1996) for IH 0707- 
495. We also included the correction for intrinsic absorption, 
0.14, the same in all sources, after Whittle (1992). The host 
galaxy luminosity of IRAS 13224-3809 was estimated on the 
basis of V band magnitude of the active nucleus 15.2 from 
Young et al. (1999) and large aperture value oi EWH/3 equal 
to 23 A given by Boiler et al. (1993). 

We see that SI galaxies and 3C 273 follow the bulge 
mass-black hole mass relation determined by Laor (1998) 
and Wandel (1999) but they are clearly shifted from the 
Magorrian et al. (1998) relation for low mass objects. For- 
mally, in our data the relation holds 

(2) 



bulge OC MbH , 



if the black hole mass is determined from the variability 
(Sect. |). 

The relation recently found by McLure & Dunlop is 
close to the original relation of Magorrian et al. (1998) and 
it is clearly inconsistent with our results. 

NLSl galaxies are strongly shifted from the SI distri- 
bution towards low black hole mass range for approximately 
the same range of bulge masses of host galaxies. 



5 DETERMINATION OF THE L/Ledd RATIO 

OF AGN 



We determine the bolometric luminosity of AGN as a sum 
of two independent contributions from X-ray and opti- 
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Table 1. Properties of AGN. Black hole masses are determined from variability (Section H), photoionization and reverberation method 
(Wandel, Peterson & Malkan 1999), accretion disk fits (Section p.2| ), and the first value is used to calculate the first luminosity to the 
Eddington luminosity ratio and the mass value from the disk fits are used to calculate the second ratio. Last column gives the mass of 
the bulge. Units: M is in Mq and L is in erg s~^. 



object 


type 


logA/— 


logMf^ 


logA/— 


logMflJ" 


logitoi 


T 1 T var 


T 1 T disk 


\0gMbulge 


NGC 5548 


81. 5 


7.39 






7.31 


44.36 


0.074 


0.089 


11.47 


NGC 7469 


81. 2 


7.03 


6.87 


6.88 


7.41 


44.30 


0.15 


0.062 


11.03 


MCG -6-30-15 


81.2 


5.94 






6.98 


43.85 


0.646 


0.059 


10.23 


NGC 4151 


81.5 


7.57 


7.35 


7.08 


7.03 
7.33 


43.70 
43.93 


0.010 
0.018 


0.037 
0.032 


11.05 


NGC 5506 


81.9 


6.86 






7.23 


43.88 


0.083 


0.036 


10.78 


ESQ 103-G35 


81.2 


7.14 






7.44 


43.39 


0.014 


0.007 


10.62 


NGC 3516 


81.5 


7.25 




7.30 


7.49 


44.46 


0.13 


0.074 


11.47 


NGC 4051 


NLSl 


5.39 


5.37 


6.15 


6.50 


42.40 


0.082 


0.006 


10.91 


IRAS 13224-3809 


NLSl 


4.52 






8.01 


44.95 


215 


0.069 


12.16 


IH 0707-495 


NLSl 


5.61 






6.77 
6.85 


43.69 
44.60 


0.95 
7.73 


0.066 
0.45 


11.39 


PC 1244+026 


NLSl 


5.85 






6.62 


44.13 


1.51 


0.26 


11.42 


PHL 1092 


NLSl 


6.09 






8.26 


45.70 


32 


0.22 


— 


3C 273 


QSO 


8.12 






8.61 


46.58 


2.3 


0.74 


12.39 


PKS 2155-304 


QSO 


7.39 






8.12 


46.21 


5.26 


0.98 


— 



cal/UV band on the basis of the integrated 2-10 keV flux 
and vFi, luminosity at 2500 A . We determine the bolo- 
metric luminosity at the basis of the following formula: 
Lj^j = 4.6xz^L^(2500A)-|-4.1xL(2-10keV). Adopted bolo- 
metric corrections were determined at the basis of the broad 
band spectrum of a Seyfert galaxy NGC 5548 (Magdziarz 
et al. 1998), although we did not include the entire unob- 
served XUV range, and of a Narrow Line Seyfert 1 galaxy 
PG 121H-143 (Janiuk, Czerny & Madejski 2001). This ap- 
proach seems to be more appropriate than using only 2-10 
keV flux with a bolometric correction 27.2 for all sources 
(Padovani & Rafanelli 1988, Hayashida et al. 1998) or us- 
ing only the optical flux (Lboi = 9AI/A (5100 A), Kaspi et 
al. 2000) since Seyfert galaxies clearly show a large range of 
the Big Blue Bump to the hard X-ray power law luminosity 
ratios (Walter & Fink 1993), particularly if NLSl galaxies 
are included in the sample. 

The data for NGC 5548 was taken from Magdziarz et 
al. (1998), X-ray 2-10 keV fluxes for MCG -6-15-30, NGC 
4051, NGC 4151, ESQ 103-G35, NGC 5506, PKS 2155- 
304 and Cyg X-3 are from Hayashida et al. (1998), for 
NGC 3516 from Nandra et al. (1997), and for IRAS 13224- 
3809, IH 0707-495, PG 1244+226 and PHL 1092 are from 
Leighly (1999a). UV fluxes for MCG -6-15-30, and 3C 273 
are from Walter & Fink (1993), NGC 4051 is from Mathur 
et al. (2000), NGC 7469 from Nandra et al. (1998). For 
PG 1244+026 the optical flux was taken from Miller et 
al. (1992), and for IH 0707-495 from Giannuzzo & Stirpe 
(1996). UV fluxes for other sources were obtained from their 
V or B magnitudes given in NED. For IH 0707 we also used 
a second set of data coming from much older observations 
when the source seemed to be fainter: optical flux and 2-10 
keV flux were taken from Remillard et al. (1986). Second 
set of data for NGC 4151 comes from the later observations 
in 1993: B magnitude is taken from Lyuty & Doroshenko 
(1999) and X-ray flux from Weaver et al. (1994). 



6 DISCUSSION 

6.1 Comparison with black hole mass 

determination from Broad Line Region 

The applicability of the method based on the Power Density 
Spectrum can be tested using those objects which black hole 
mass was determined independently. 

Some of the sources in our study were extensively moni- 
tored (NGC 4151, NGC 4051, and NGC 7469) which allowed 
to use reverberation or photoionization method, or both, to 
determine the mass of the central black hole. The values 
obtained by Wandel, Peterson & Malkan (1999) are given 
in Table [D. NGC 3516 was also monitored and the black 
hole mass determined by reverberation method was given 
by Wanders & Home (1994). 

We see that for three objects there is a good agreement 
between the black hole mass determined from X-ray vari- 
ability with the values from both the ionization and rever- 
beration method. Methods based on the Broad Line Region 
properties are themselves accurate within a factor of three, 
according to the discussion recently made by Krolik (2000) , 
and the mass from PDS is contained in the expected range. 
In case of NGC 4151 the black hole mass determined from 
variability is larger by a factor of 1.7 than the value from 
photoionization and by a factor of 3 than the value from re- 
verberation. In that case the power density spectrum taken 
from Hayashida et al. (1998) may not be determined ac- 
curately. This object, having larger mass, requires longer 
monitoring than NGC 4051 and the Ginga results may not 
be adequate. The source was earlier observed by EXOSAT 
and the power spectrum was determined by Papadakis & 
McHardy (1995) but without the normalization. If we nor- 
malize their spectrum according to the provided mean in- 
tensity and the excess variance and use their PDS we obtain 
even larger black hole mass (log Mbh = 8.43). However, 
true high frequency part of the PDS may be actually flatter 
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Figure 4. The relation between the mass of the central black 
hole Mg°^ determined from variability and the mass of the host 
galaxy bulge: filled squares show SI galaxies and 3C 273, and 
open pentagons mark NLSl. Continuous line shows the relation 
of Wandel (1999) for Seyfert galaxies, dashed line is for quasars 
and normal galaxies (Laor 1998), and dot-short-dash line is the 
original relation of Magorrian et al. (1998). Star marks the posi- 
tion of the Milky Way. 



than the slope -2 given in the paper due to the power leakage 
effect discussed by the authors. If we arbitrarily assume that 
the turning point is correct but the slope is -1.5 we obtain 
log Mbh = 7.5. This last value is closer to the values deter- 
mined from the BLR properties, but this discussion mostly 
points out a need of high quality PDS of AGN. 

6.2 Comparison with black hole mass 

determination from accretion disk fitting 

For most of the sources from our sample such determina- 
tions are not available (and for other objects from Wandel 
et al. (1999) sample the power spectra were not determined). 
Therefore, in order to test better the method we determine 
the black hole masses for all sources based on the following 
assumptions: (1) outer part of the flow responsible for the 
UV emission is well modeled by the standard accretion disk 
(2) the X-ray emission is produced closer in, and the total 
efficiency of conversion of the accreting mass into observed 
radiation is 1/16, as appropriate for non-advecting flow onto 
a Schwarzschild black hole. 

The first assumption results in adopting the relation 
between the black hole mass, accretion rate, intrinsic lumi- 
nosity Li, at a frequency v after Tripp, Bechtold & Green 
(1994) 



logM + logM = 1.5{logL„ ~ 19.222) - 0.5logv, 



(3) 



Figure 5. The relation between the mass of the central black 
hole Mg°^ determined from variability and the mass determined 
on the basis of accretion disk theory MgS!': filled squares show 
SI galaxies and QSOs, and open pentagons mark NLSl. Dotted 
line marks the expected relation M'^^ 



Mflj'. 



This relation, combined with the second assumption 
and calculated at 2500 A, gives 

logM = 1.5log{uLf°°) - logLtoi ~ 13.44. (4) 

Since it relies on disk assumption we call it Mg^j' in Table III 
Such an approach does not take into account the depar- 
ture of the locally emitted spectrum from the black body, 
and advection or outflow may lead to smaller energetic ef- 
ficiency 77 than assumed. The resulting black hole mass de- 
pends on those factors as 



M xfr] 



(5) 



(Janiuk et al. 2001), where / is the color temperature to 
the effective temperature ratio so the mass may be either 
underestimated or overestimated by a factor of a few. 

The relation between the masses determined in two 
ways are shown in Fig. H. 

The black hole masses obtained in this way for SI galax- 
ies and QSOs are quite similar to values obtained from vari- 
ability. Mean value of the logarithm of the black hole mass 
for a SI galaxy is equal to 7.03 from variability and 7.27 
from disk fitting, and for the two quasars these two numbers 
are equal 7.75 and 8.36, correspondingly. It might indicate a 
slight systematic non-linearity required in the formula con- 
necting the mass with the frequency t'o.ooi (see Section 2.1) 



where M is in solar masses, M is solar masses per year, f 
in Hz, and L^ in erg s~^ Hz~^. 



but the uncertainties are currently too large to address this 
issue. 

Recent paper of Collin & Hure (2001) claimed a con- 
siderable problem with accretion disk luminosity but we did 
not found any evidence for such an effect (see Appendix A). 

However, we note a large discrepancy between the meth- 
ods in the case of NLSl, with mean masses equal 5.49 from 
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variability and 7.16 from disk fitting. The power spectra for 
four out of five objects comes from Hayashida (2000) and 
they are of low quality, as indicated by large errors given 
by the author, but better spectra are not available. Higher 
quality spectrum available for NGC 4051 gives the value of 
the mass located between the two values provided by the 
BLR studies which may indicate that the scaling with NGC 
5548 which is SI may be nevertheless appropriate for NLSl 
objects. In this case a larger error may be actually connected 
with the way how we determined the bolometric luminosity 
which influenced both the luminosity to the Eddington lumi- 
nosity ratio as the determination of the Mg^'' . Our formula 
(given in Sect, a) is appropriate for sources either strong 
in X-rays (as SI galaxies) or sources with a Big Blue Bump 
having a maximum in vF^, in UV band, like PG 1211-1-143 or 
quasar composite spectrum of Laor et al. (1997). However, 
for some NLSl the Big Blue Bump may be located in very 
far UV or even soft X-rays. One example of such a source 
(RE J1034-I-396) was found by Puchnarewicz et al. (1995). 
In this extreme case the bolometric luminosity estimated 
from our formula is 6.4 x 10*'^ erg s~^, while direct integra- 
tion of the broad band spectra would give about 2 x 10^^ 
erg s~^, a factor of 30 higher, and the first value would give 
a black hole mass by a factor of 30 too high. Large error is 
therefore possible, particularly for sources which are weakly 
variable in the optical band like NGC 4051, which may also 
have a Big Blue Bump shifted considerably towards higher 
energies, although not as high as in RE J1034-I-396. There- 
fore, Mg^ niay be sometimes too high while on the other 
hand the values of Ml% for IRAS 13224-3809, IH 0707- 
495 and PHL 1092 are clearly too low, leading to extreme 
super-Eddington luminosities. 

Recent results for the X-ray variability of NLSl galaxy 
Akn 564 does not solve the issue since only the rising and the 
flat part of the power spectrum on :/ x Power diagram were 
detected (Pounds et al. 2001) thus giving only a lower limit 
on the frequency t'o.ooi (~ 10~* Hz and subsequently an 
upper limit on the black hole mass (IO^Mq) since the mon- 
itoring procedure did not allow to study higher frequencies. 

Summarizing, the PDS method seems to lead to good 
black hole mass estimation for normal Seyfert galaxies and 
QSOs the high quality PDS is available. The conclusion is 
less firm for NLSl but this is probably caused by low quality 
of the available PDS spectra for those objects. 

6.3 Luminosity to the Eddington luminosity ratio 

The luminosity to the Eddington luminosity ratio for SI 
galaxies is ~ 0.05 up to 0.07, practically independently from 
the method used to determine the masses of their central 
black holes. 

Only one source, MCG -6-30-15, has this ratio higher 
than 0.15. This object is rather exceptional with respect to 
the properties of its X-ray spectra: it has relatively broad 
iron Ka line, with a strong red wing (Tanaka et al. 1995, 
Nandra et al. 1997). MCG -6-30-15 has also very soft spec- 
trum, with a hard X-ray photon index in ASCA data 1.95 
and 2.04 in two data sets (Nandra et al. 1997). The Ha line 
of MCG -6-15-30 is rather narrow, about 2000 km s~^, and 
the Hp line also seems quite narrow (Morris & Ward 1988), 
although the exact value may strongly depend on the de- 
composition of the line into narrow and broad component. 



Turner et al. (1999) quote the value of 1700 km s~^ after 
Pineda et al. (1980). Classification of this source should be 
perhaps reconsidered. 

NLSl galaxies show systematically higher luminosity 
to the Eddington luminosity ratios (although the strength 
of this effect depends on the black hole mass measurement 
method). This is consistent with suggestion of Pounds, Done 
& Osborne (1995), confirmed by a number of later studies 
(e.g. Czerny, Witt & Zycki 1996, Wandel 1997, Brandt & 
Boiler 1998, Kuraszkiewicz et al. 2000). 

Black hole mass estimate based on the variability 
method leads for three sources to unrealistic, highly super- 
Eddington values of the accretion rate. Higher quality power 
spectra for these sources are clearly needed. Spectral fitting 
approach lead to much smaller ratios, with mean value 0.18, 
more than three times higher than for SI galaxies. NLSl 
galaxies have on average higher luminosity to the Eddington 
luminosity ratio although some NLSl galaxies have those ra- 
tios (usually only slightly) smaller than some of SI galaxies. 

Only NGC 4051 among NLSl has an exceptionally low 
luminosity to the Eddington luminosity ratio. The mass de- 
termined from the variability seems to be correct for this 
source, as it is intermediate between the other values quoted 
in the literature: Hp study gave (1.1 - 1.4) x IO'^Mq, (Pe- 
terson et al. 2000), Kaspi et al. (2000) obtained (1.3-1.4) x 
10® M©, and Wandel et al. (1999) values were 2.3 x 10^ Mq 
and 1.4 X IQ^Mq, as given in Table 0. It may be that the 
source bolometric luminosity is by an order of magnitude 
underestimated if the broad band spectrum is dominated 
by the unobserved XUV range. On the other hand Peterson 
et al. (2000) argue that this source may have relatively nar- 
row Hf} line due to high inclination (~ 50°, Christopoulou 
et al. 1997) instead of high L/LEdd- Better determination of 
the PDS for this source may resolve this issue. 

Two QSO included in our sample have rather high 
L/LEdd ratios. The black hole mass of 3C 273 based on 
variability and on our simple spectral fit is smaller than the 
value determined from detailed spectral fitting (e.g. Leach, 
McHardy & Papadakis 1995). The problem may be con- 
nected with too short monitoring, i.e. the lack of coverage 
of frequencies close to flattening point and subsequently, a 
problem with normalization of the PDS. On the other hand, 
in the Leach et al. (1995) flts the energy required to heat the 
Comptonizing corona is not included in the computation of 
the accretion rate, possibly giving too low accretion rate and 
too high mass. The analysis of longer timescale variability 
in these sources is clearly needed. 



6.4 Spectral states of AGN 

There are considerable similarities between the accretion 
process onto stellar black holes and supermassive black 
holes: similar range of X-ray spectral slopes (Zdziarski 
1999), similar correlation between the normalization of the 
reflected component and X-ray spectral slope (Zdziarski, 
Lubihski & Smith 1999) and basically similar featureless 
PDS (Hayashida et al. 1998). The shape of X-ray PDS for 
both NLSl galaxies and normal SI galaxies look quite sim- 
ilar to the GBH in their hard state. 

Detailed comparison of AGN and GBH shows, how- 
ever, a systematic difference: (i) GBHs in their soft state 
show very little variability, with rms at a level of a few per 
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Figure 6. The comparison of the optical power density spectrum 
times frequency for NGC 5548 (Czerny et al. 1999; disconnected 
points) with the X-ray power density spectrum of NGC 5548 from 
Chiang et al. (2000) (histogram: continuous line - XTE; stars con- 
nected with continuous line - fit) and for galactic black holes: Cyg 
X-1 in hard state from Revnivtsev et al. (2000) (open squares 
connected with dot line), Cyg X-1 in soft state from Gilfanov et 
al. (solid triangles connected with long dash line). Power den- 
sity spectra of NGC 5548 were shifted by 6.83 in logiy) which 
corresponds to rescaling them to a black hole mass of WMq . 



cent (van der Klis 1995) while even disc-dominated AGN 
like quasars generally show higher variability both in opti- 
cal (Giveon et al. 1999) and X-ray band (e.g. Yaqoob et al. 
(1994) for PG 1211). 

This apparent inconsistency can be qualitatively under- 
stood in the following way. The similarity of the X-ray PDS 
for all kinds of AGN reflects the fact that this PDS is usually 
calculated using the 2-10 keV energy band which is always 
dominated by the hard X-ray power law. X-ray emission of 
the galactic sources in their soft state is, however, domi- 
nated by the soft component which roughly corresponds to 
the optical/UV component in AGN. 

The same variability mechanism responsible for the X- 
ray PDS of AGN and GBH in their hard states does not 
mean that the two must be identical. Actually, we would 
rather expect slight systematic shift between AGN and GBH 
(apart from the one resulting from the mass ratio) since the 
observed PDS of galactic sources show some dependence of 
PDS on the studied spectral range (e.g. Cui et al. 1997, 
Nowak et al. 1999, Gilfanov, Churazov & Revnivtsev 1999). 
The spectrum of NGC 5548 is also dominated by hard X- 
ray power law, without a strong contribution from soft X- 
ray excess. However, soft seed photons for Comptonization 
have much lower energy in case of AGN and have to undergo 
more scatterings before reaching the same energy as photons 
in Cyg X-1. Therefore, even if the hot medium temperature 
and geometry is identical in an AGN and a GBH the Comp- 



tonization process is not identical and an AGN PDS may 
resemble more the PDS of a GBH at somewhat higher en- 
ergy band. Numerical simulations should be used to study 
such a scaling properties with the black hole mass, i.e. with 
the soft photon temperature. 

The essential differences between AGN and GBH are 
reflected in the behavior of their soft spectral components. 
Cyg X-1 in disk-dominated soft state shows very moderate 
variability. The behavior of disk component in NGC 5548 is 
well represented by the optical PDS shown in Figure g. The 
errors in determination of the optical PDS are large but we 
see that the overall shape is more similar to that of Cyg X-1 
in its hard state than the X-ray PDS. It may mean that the 
basic variability mechanism is similar in both objects and 
that the spectral state of NGC 5548 really corresponds to 
low/hard state, as suggested by its L/ Lsdd ratio, although 
the normalization of the optical PDS of NGC 5548 is much 
higher than that of Cyg X-1. 

Such an effect is not surprising. According to a standard 
Shakura-Sunyaev disc model (Shakura & Sunyaev 1973) we 
deal with the different ratio of the radiation to the gas pres- 
sure in accretion discs in AGN and GBH at the same L/LEdd 
ratio. Time evolution and hot plasma formation rate can 
clearly depend on this ratio. Rozariska & Czerny (2000) 
show the expected difference between AGN and GBH accre- 
tion flow stability in the case of accreting corona model. In 
this case an evaporation branch exists which limits optically 
thick solutions to high accretion rates only, with a transition 
to an advection-dominated accretion flow (i.e. ADAF; Ichi- 
maru 1977, Narayan & Yi 1994, for a review see Narayan, 
Mahadevan & Quataert 1998) for low accretion rates. The 
accretion rate at which the evaporation branch forms does 
not depend on the mass of the black hole but the extension 
of the gas pressure dominated branch of a standard disc does 
depend on the mass considerably. As a result, AGN never 
achieve a stable high state since before the accretion rate 
drops enough to allow the gas pressure to dominate the disc 
evaporates. 

6.5 Black hole to galactic bulge mass ratio 

The mass of the central black hole residing in local mas- 
sive galaxies is well correlated with the mass of the bulge 
(Magorrian et al. 1998). It suggests a close relation between 
the host galaxy and the formation, or evolution, of its central 
black hole. Laor (1998) showed that nearby quasars studied 
by Hubble Space Telescope follow the same relation. 

However, Wandel (1999) studied this relation for Seyfert 
galaxies and found that they form a similar relation as nor- 
mal galaxies and quasars but systematically displaced to- 
wards small black hole masses. Similar study performed by 
Mathur et al. (2000) for Narrow Line Seyfert 1 galaxies 
showed even further displacement towards small black hole 
masses. 

Our mass determination used in Fig. H was based on 
PDS analysis so it was different from method used by 
Mathur et al. (2000) and Wandel (1999), and mostly per- 
formed for other objects. Our results, however, confirm the 
trends found by those authors - Seyfert galaxies follow basi- 
cally the relation determined by Wandel but NLSl galaxies 
depart from it systematically. The strength of the effect, 
however, significantly depends on the adopted method of 
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ing periodically temporary enhancement or suppression of 
the accretion it leads to transient behavior of many sources 
in timescales of years (e.g. King 1995). It was suggested 
that similar instability may operate in AGN (e.g. Clarke 
& Shields 1989, Siemiginowska, Czerny & Kostyunin 1996, 
Hatziminaoglou, Siemiginowska & Elvis 2000) although ex- 
pected timescales are of thousands to millions of years. Such 
instability may modify the accretion rate in the central parts 
of the flow thus being responsible for temporary change 
of the object status from highly active (NLSl) to regular 
(Seyfert 1) and perhaps quiescence (LINERS) in similar 
timescales. On the other hand it is not clear that this in- 
stability operates in AGN since the self-gravity effects are 
extremely important in the outer partially ionized parts of 
discs around massive black holes, leading to significant mod- 
ification of the disc structure in comparison with discs in low 
mass objects (e.g. Hure 1998). The fact that NLSl galaxies 
and normal Seyfert 1 galaxies have different black hole mass 
to bulge mass ratio and therefore different evolutionary sta- 
tus argues against the temporary transition between NLSl, 
SI and LINER stages and in favor of the self-gravity pre- 
venting the ionization instability to operate in AGN. Further 
studies, however, are needed to confirm this conclusion. 



Figure 7. The relation between the mass of the central black hole 
determined from spectral estimates {Mgfj' from Table hi) and the 
mass of the host galaxy bulge: filled squares show SI galaxies and 
3C 273, and open pentagons mark NLSl. Continuous line shows 
the relation of Wandel (1999) for Seyfert galaxies, dashed line is 
for quasars and normal galaxies (Laor 1998), dot-short-dash line 
is the original relation of Magorrian et al. (1998), and long dash 
line is of Laor (2001). Dotted line shows our relation based of 
equation (H). Star marks the position of the Milky Way. 



black hole mass measurement for NLSl galaxies: NLSl are 
shifted by ~ 1.5 orders of magnitude in Fig. Ml (black hole 
mass from variability) and it is marginally shifted by ~ 0.5 
magnitude in Fig. M (black hole mass from spectral fits). 
We therefore may even neglect the difference between NLSl 
galaxies and SI galaxies and look for a universal relation be- 
tween the bulge mass and black hole mass for our objects. 
Spectral black hole mass determination gives 



M, 



bulge 



oc AfR 



(6) 



closer to the non-linear relation found recently by Laor 
(2001). However, without NLSl galaxies this relation re- 
mains linear for our objects even in the case of the black 
hole mass determination from spectral fits. 

The quasar 3C 273 in our analysis follows the relation 
for Seyfert galaxies. 

Recent study of normal galaxies seem to indicate that 
even for these objects the black hole mass and bulge mass 
relation is not universal, showing trends with the mass of 
the galaxy (Merritt & Ferrarese 2000). 



7 CONCLUSIONS 

• high frequency part of the PDS is luminosity- 
independent and can be used to determine the mass of the 
central object 

• Cyg X-1 may be used as a reference object for AGN 
although systematic differences between AGN and GBH in 
the appropriate normalization of the PDS up to an order of 
0.4 in logarithm are not excluded 

• NLSl galaxies may possibly require another scaling 
than SI galaxies 

• the L/LEdd ratio based on mass determination from 
PDS has large values for NLSl galaxies and lower values for 
most BLSl galaxies. The effect is still present, although sig- 
nificantly reduced, if complementary black hole mass mea- 
surement based on simple spectral fitting is adopted, 

• NLSl galaxies seems to have systematically lower black 
hole masses in comparison to bulge masses of the host galax- 
ies which supports the conclusion of Mathur (2000) about 
their different evolutionary status. MCG -6-30-15 occupies 
its appropriate part of the diagram despite having relatively 
large L/LEdd, 

• ionization instability operating in cataclysmic variables 
and X-ray novae probably cannot be responsible for varia- 
tions of accretion rate corresponding to transitions between 
NLSl and BLSl state over the periods of thousands - mil- 
lions of years. This instability may be suppressed in AGN 
by self- gravity effect. 



6.6 Ionization instability in AGN 

Ionization instability, first found to be responsible for the 
outbursts of cataclysmic variables (e.g. Meyer & Meyer- 
Hoffmeister 1984, Smak 1984) operates also in binary sys- 
tems containing accreting neutron stars or black holes. It 
develops in outer, partially ionized part of the disc. Caus- 
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APPENDIX 

In their recent paper Collin & Hure (2001) analyzed the 
objects from the sample of Kaspi et al. (2000) and they found 
serious deficiency of the putative accretion disk luminosity 
in the optical band. Therefore, we checked our approach 



outlined in Section 6.2 against this sample. 



If we assumed the same prescription for the bolometric 
luminosity as Kapis et al. (2000), i.e. 

Ltoi = 9iyL,{5100A), (7) 

determine the accretion rate at the basis of the accretion 
efficiency rj 

2 



Lboi = vMc 



(8) 



and adopt the value of the black hole mass from either of 
the two fits (mean or rms. Equation 9 and 10 of Kaspi et al. 
2000), we can calculate the expected monochromatic accre- 
tion disk luminosity at 5100 A, in order to compare it with 
the observed value. 

Expressing M in solar masses and M in solar masses 
per year we obtain from Tripp et al. (1994) (see Equation ^ 
in case of mean mass fit 



1.031og!/L°'"'(5100)- 



1.22--log(16r,),(9) 



log//L*'''=(5100) 

and for rms fit 

logi/Lf'''=(5100) = 0.941ogi/Lf ''(5100)-h2.56-^log(167?),(10) 



and for sources at the typical luminosities i/Li, 



5100 



10^^ erg 



s~^ there is no discrepancy between the disk luminosity and 
the actual luminosity, contrary to Collin & Hure (2001), if 
r] = 1/16 and cosi = 0.5 is assumed, as we did in our consid- 
eration. Collin & Hure (2001) adopted 77 = 0.1 and cosi = 1 
which partially although not fully explains their conclusion 
about too low disk luminosity. Applying the top view may 
actually be correct; however in that case such a factor should 
be probably included in the formula for the isotropic bolo- 
metric luminosity as calculated from the observed flux which 
will weaken the dependence to the power 1/3. 
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